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A low dimensional tool for flow-structure interaction problems based on Proper Or- 
thogonal Decomposition (POD) and modified Linear Stochastic Estimation (mLSE) has 
been proposed and was applied to a Micro Air Vehicle ( MAV ) wing. The method uti- 
lizes the dynamic strain measurements from the wing to estimate the POD expansion 
coefficients from which an estimation of the velocity in the wake can be obtained. 

For this experiment the MAV wing was set at five different angles of attack, from 
0° to 20°. The tunnel velocities varied from 44 to 58/f/sec with corresponding Reynolds 
numbers of 46,000 to 70,000. A stereo Particle Image Velocimetry ( PIV ) system was used 
to measure the wake of the MAV wing simultaneously with the signals from the twelve 
dynamic strain gauges mounted on the wing. 

With 20 out of 2400 POD modes, a reasonable estimation of the flow field was ob- 
served. By increasing the number of POD modes, a better estimation of the flow field 
will occur. Utilizing the simultaneously sampled strain gauges and flow field measure- 
ments in conjunction with mLSE , an estimation of the flow field with lower energy modes 
is reasonable. With these results, the methodology for estimating the wake flow field 
from just dynamic strain gauges is validated. 


Introduction 

The Micro Air Vehicle (MAV) program was started 
from a Defense Advanced Research Projects Agency 
(DARPA) initiative. DARPA wanted an aircraft with 
a wingspan of no more than six inches to be developed. 
Government, industry and academic institutions have 
worked on the problem and each has their own purpose 
for development. MAVs, in a military role, can be used 
to show foot soldiers what is over the next hill. For 
surveillance purposes, the MAV can be used to observe 
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without, being seen. For a scientific application, the 
MAV can be used for the cost-effective remote sensing 
of atmospheric conditions . 1 

One problem with all MAVs is the low Reynolds 
number flow that they experience while in flight. Since 
most MAVs fly in atmospheric conditions, gust and 
turbulence can affect the flight and controllability of 
these vehicles. By understanding how the MAVs will 
react to gust and turbulence, control systems can be 
developed to reduce the effects of the external load- 
ings. Developing a control system for the motion of the 
MAV’s wing and its wake flow field requires them to be 
simultaneously sampled because the problems are cou- 
pled. With a feedback control system, the enhanced 
MAV wing might have improved flight characteristics, 
such as increased maneuverability or stability. In the 
work presented here the focus is on the sensing aspect 
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of the control problem. We use surface strain mea- 
surements to estimate the POD expansion coefficients 
in the wake, which are then used to estimate the ve- 
locity wake flow field when combined with the POD 
eigenfunctions . 

Background 

In 1984, Sneyd 2 derived the aerodynamic coefficients 
and stability of a sail with no slack while it was at- 
tached to a rigid structure, i.e. fixed leading and 
trailing edges where camber is accommodated by the 
stretching of the fabric. He found that sails with no 
slack have the important property of neutral static 
stability and if the wings are sufficiently stretchable, 
they may have a slight degree of positive stability. 
This property must be important in understanding the 
flight dynamics of bats and pterodactyls whose wings 
consist of a tensioned elastic membrane without slack. 

In 1995, Smith and Shyy 3 performed computations 
of the harmonically forced, unsteady viscous flow over 
a flexible, two-dimensional membrane wing. The sim- 
ulation of a harmonically forced free-stream flow also 
proves to be a useful vehicle for demonstrating some of 
the more generic features of membrane wing mechan- 
ics. 

Smith and Wright 4 used computations of the har- 
monically forced, unsteady incompressible flow over 
an extensible flexible membrane airfoil. The role of 
viscosity was shown to be significant in the harmoni- 
cally forced unsteady flow about a membrane wing. In 
the unsteady scenario, the influence on viscosity is en- 
hanced since the acceleration and deceleration of the 
free stream velocity strongly influences the separation 
and reatt.achment of the flow. 

In 2001, Dowell and Hall 5 showed various physical 
models for a fluid undergoing time-dependent, motion, 
the distinction between linear and nonlinear models, 
time-linearized models and their solutions in the fre- 
quency or time domain. They extended treatment 
of the modal character of time-dependent, flows and 
the construction of reduced-order models base on an 
expansion in terms of fluid modes. Other works by 
Naudascher and Rockwell 6 and Dowell 7 on this sub- 
ject offer better understanding of problems that, deal 
with the flow and structure. 

Lumley, 89 developed the Proper Orthogonal De- 
composition (POD) technique as a tool to identify 
the most energetic coherent, structures, or eddies , con- 
tained within a turbulent, flow field. The first, valida- 
tion of the POD was performed by two of his graduate 
students. Payne 10 found that, the energy associated 
with the “dominant, eddy” in the wake of a cylinder 
was not significantly greater than any of the other 
structures. Bakewell and Lumley 11 then applied the 
POD to the near wall flow region of a pipe flow and 
showed that, the “dominant, eddy” contained roughly 
90% of the energy. These initial efforts suggested the 


power of the POD could be a useful tool for the study 
of turbulent, flows. 

In 1998 Bonnet, et al .' 12 ,9 provided a comparison of a 
few tools for the decomposition and analysis of turbu- 
lent flows that, placed an emphasis on the detection of 
large scale coherent, structures. Some techniques that, 
were covered are Proper Orthogonal Decomposition, 
Linear Stochastic Estimation (LSE) and Conditional 
Sampling Approaches. 

In 2002, Taylor and Glauser 13 performed an experi- 
ment with a separating shear flow in the ActiveWING 
facility at Clarkson University. This facility is a two- 
dimensional backward facing ramp with an adjustable 
flap above the ramp that, creates conditions from fully 
separated to fully attached flow on the ramp. Parti- 
cle Image Velocimet.ry (PIV) and wall pressure were 
simultaneously sampled at. a number of Reynold num- 
bers and many angles of attack. Moreover, the instan- 
taneous flow field can be reconstructed for any flow 
condition measured using the instantaneous pressure 
measurements on the wall. 

MAV Wing 

The dimensions of the MAV wing used for the ex- 
periment are shown in Figure 1. The span and chord 
lengths of the wing are 5.75 and 3 inches, respectively. 
The aluminum frame consisted of one spar and four 
ribs. The frame has a 0 0 camber. An inch of the 
frame was used to hold the wing into the mounting 
block. The mounting block was then attached to the 
tunnel floor. The wing’s angle of attack can be ad- 
justed from —20° to 20° by rotating the mounting 
block. To reduce the influence of the tunnel boundary 
layer, the wing was mounted 1/2 inch above the tunnel 
floor. 



Fig. 1 Dimensions of MAV Wing 


A rectangular latex membrane 3 by 5.25 inches, was 
used for the skin of the wing. Three membrane regions, 
2.625 by 1.25 inches, were formed between each pair of 
ribs. Super 77™ a spray adhesive from 3M™ was 
used to adhere the latex membrane to the frame. 
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Experiment 

The experiment was conducted in the NASA Lan- 
gley 15-Inch Low Speed Wind Tunnel. This tunnel is 
a low-speed, closed return, atmospheric facility used 
primarily for fundamental flow physics research. The 
maximum flow speed in the test section is approxi- 
mately 115 ft/sec. All data included here was taken in 
the test section closest to the inlet. 

A Laranko Stereo Particle Image Velocimetry, 15 
PIV, was employed to capture the wake flow field of 
the MAV wing. The system uses a 400m J YAG Laser 
to illuminate the seed particles in the flow field. Two 
2Kx2K digital cameras are used to capture the flow. 
Mineral oil was atomized and was used for the seed 
particulate which has the diameter of .1 to 5 /im’s and 
can follow the flow up to the frequency of 1 kHz. 

Using I-Deas Finite Element Modeling Program the 
placement of the strain gauges were determined using a 
Normal Mode Dynamic Solver. A total of six locations 
of high strain have been determined. At each of these 
locations a bi-axial strain gauge were placed. Figure 2 
shows the placement of all twelve strain gauges, they 
are numbered for identification purposes. 

Placement of Strain Gauges for MAV wing Drawing 01 



Fig. 2 Placement of Strain Gauges on the Frame 

Each strain gauge used has a gauge length and width 
of 3mm ' s and of 2 mm's, respectively. The matrix 
length and width are 7mm s and 4mm' s, respectively. 
The Gauge Factor of 2. 120 ±0.5% at 24° ( 7 . The resis- 
tance is 350.0 ± 0.3% ohms at 24°C with a Transverse 
sensitivity of (±0.7 ± 0.2)% at 24°C. There are three 
wire tabs so that temperature compensation can be 
taken into account. The strain gauges were placed on 
the aluminum with an adhesive. The wires are then 
lightly adhered to the wings so that they will not vi- 
brate in the flow. 

To determine the instantaneous wake flow field of 
the MAV wing, the PIV system and strain gauges 
were used. Because the PIV system is independent of 
the Zonic data acquisition system, used for the strain 
gauge measurements, a method to connect the systems 
had to be developed, enabling the simultaneous sam- 
pling of the strain gauges and the wake flow field to 


occur. The methodology that, was developed involved 
t aking the Q-swit.ch signal from the first. PIV laser and 
allowing the DAQ system to sample the signal and can 
be seen in Figure 3. The Q-swit.ch signal works on a 
Transist.or-Transist.or Logic (TTL) pulse of ±5 volts 
10 Hz square wave, that, is sent, from the control box 
to the lasers. It. can also be started with an external 
trigger that, is set. with a TTL pulse from a signal gen- 
erator. Every time a ±5 volt, pulse is sent, from the 
controller to the lasers, the lasers will fire. A trigger 
setting of a negative slope and 50% change in the full 
scale voltage was used t.o start, the sampling of the 
strain gauges. Because the Q-swit.ch is being sampled 
simultaneously with the strain gauges, the PIV images 
are then considered to be simultaneously sampled with 
the strain gauges. 


Signal Generator 



Q-Switch 2 

Fig. 3 Wire Diagram of Simultaneous Sampling 
Setup 

For the experiment the MAV wing’s angle of at- 
tack was changed from 0° to 20° at 5° intervals. The 
tunnel’s air speed was changed from 44 ft/s to 58 ft/s 
which gives a Reynolds number of 40, 000 and 70, 000 
respectively. The PIV images were captured at c/2 
and 2c/3, 1.5 and 2.0 inches, respectively downstream 
of the wing. With five different angles of attack, two 
tunnel airspeeds, and two PIV image planes, a total of 
20 runs are required for this experiment. 

Analysis 

In 1967 Lumley 8 proposed the Proper Orthogo- 
nal Decomposition, POD, as an unbiased technique 
for studying coherent structures in turbulent flows. 
The POD is a mathematical approach based on the 
Karhunen-Loeve expansion to decompose the fluctu- 
ating velocity field into a finite number of empirical 
functions. 

Adrian 16 proposed that stochastic estimation could 
be applied to unconditional data. This method used 
what Adrian called a “conditional eddy” to form a 
technique for estimating the flow field given the pres- 
ence of this defined “conditional eddy.’’ Cole et al. 17 
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then proposed that the instantaneous velocity field 
could be used, instead of defining an arbitrary indi- 
cator such as Adrian’s “conditional eddy.” 

In 1994, Bonnet et al. ls presented the Complemen- 
tary Technique which combines the LSE and the POD 
into a tool for the development of time varying low- 
dimensional descriptions of the flow. 

Taylor and Glauser 13 proposed the use of Global 
POD to reconstruct the instantaneous flow field from 
the instantaneous wall pressure measurements in a flow 
over a backward facing ramp with a cyclic flap over 
the ramp which allows for a dynamically changing ad- 
verse pressure gradient. In 2002, Schmit 19 extended 
the work of Taylor to reconstruct the instantaneous 
velocity wake flow field from the instantaneous motion 
of the MAV wing. For this experiment, twelve judi- 
ciously placed strain gauges and a stereo PIV system 
were used to measure simultaneously the instanta- 
neous strain and three components of velocity in the 
wake of the MAV. 

The analysis methodology can be summarized as: 

1. From a 2D — 3 Component velocity measurement, 
determine the two-point correlation tensor, Equation 
1, at each PIV stream-wise location, *o, measured in 
the wake. 


M 




Ui{xp,y, z),Uj(xo,y\ z 1 ) 

M 


(!) 

With Equation 1 and following the Hilbert-Schmidt 
theory (Lumley 8 ), a discrete set of solutions, Eiquation 
2, can be obtained if random fluctuations occur over a 
finite domain. 



y, y',z, z')<j)\ n \x o, y' ,z')dy'dz' = 


A (n) (®o)^ n) (*o.2/.2) ( 2 ) 

Equation 2 is solved numerically for the eigenfunc- 
tions of the two-point correlation. The can then 
be used to reconstruct the original fluctuation quanti- 
ties with Equation 3. 


m 

Ui(x 0 ,y, z, t) = (x 0 ,t)</)\ n) (xo,y,z) (3) 

n — 1 

The POD expansion coefficients of the flow, 
a n (xo,t) can be determined with Equation 4, where 
the (f)^* denotes the complex conjugate of . 

an(xo,t) = / Ui(x 0 ,y, z,t)<j) ( A ) *(xo,y, z)dydz (4) 
Jd 

2. Using the simultaneously sampled strain gauge 
measurements and the POD expansion coefficients, 


the LSE coefficient, A nk , can be determined for each 
streamwise location by: 
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3. Finally, estimate the POD expansion coefficients 
from the temporal stain gauge measurements, shown 
in Equation 6, and combine them with the eigenfunc- 
tions, to produce the instantaneous estimations of flow 
field as expressed in Equation 7. 

q 

a n est (x 0 J) = '^2A nk ei ! {t) (6) 

k = 1 

m 

Ui(x 0 ,y,z,t ) = a n est (x 0 , t)^ n: (x 0 , y, z) (7) 

n=l 


Results 

The stereo PIV system used for this experiment can 
measure the instantaneous wake flow field in a two- 
dimensional plane with three components of velocity. 
Figure 4 shows an instantaneous snapshot of the flow 
at 44/ps, downstream at c/2 and an angle of attack 
of 15°. The reader is looking downstream, in the x 
direction, from the MAV wing to the yz plane. Figure 
4 contains four separate sets of data which can be seen 
in Table 1. 




Fig. 4 Instantaneous Snapshot of the Wake Flow 
Field 


Eigenvalues and Eigenvectors 

Once the two-point correlation tensor has been 
formed the eigenfunction of the tensor can be solved 
by using Equation 2. One method used to understand 
the eigenvalues and their importance is to examine the 
energy contained up to N POD modes as depicted in 
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Upper Left U Component 

Upper Right V Component 

Lower Left W Component 

Lower Right V and W Components 

Table 1 Results Setup for Instantaneous Snapshot 
of the Wake Flow Field 

Equation 8. Note that in Equation 8, N will eventually 
go to m the total number of POD modes. 

X> (n) (8) 

n= 1 

We solved the eigenvalue problem with and with- 
out the mean velocity flow field. Figure 5 shows the 
summation of energy for each set of conditions with 
the mean included. Figure 6 shows the summation 
of energy for each set of conditions without the mean 
included. 

As seen in Figure 5 the mean mode has about 95 to 
97% of the total energy of the flow, this is out of 2400 
POD modes. From Figure 6 the summation of the tur- 
bulent modes show that the 1 st mode has about 20% 
of the total energy, this is out of 2399 POD modes. In 
both cases, as N increase in Equation 8 the total en- 
ergy available in that mode is less than the one before. 
It takes approximately 400 POD modes in both cases 
to obtain nearly 100% of the energy available in the 
flow. However, with only 20 POD modes, we obtain 
more than 99% of the total energy when we include the 
mean as seen in Figure 5 and more than 70% of the 
total energy when we examine the fluctuations only as 
shown in Figure 6. These results give a good indication 
that this flow has a low-dimensional description. 


Summation of Energy, TS =44 fps, Plane = c/2 Summation of Energy, TS =58 fps, Plane = c/2 



0 100 200 300 400 0 100 200 300 400 

POD Mode POD Mode 

Summation of Energy, TS =44 fps, Plane = 2c/3 Summation of Energy, TS =58 fps. Plane = 2c/3 


0 100 200 300 400 0 1 00 200 300 400 

POD Mode POD Mode 

Fig. 5 Summation of Energy per POD mode with 
mean included. 

Figures 7 to 10 show the first 4 POD eigenvectors, 

. with the mean included. Because the 1 st POD 
mode contains up to 95 % of the total energy it should 


Summation of Energy, TS =44 fps, Plane = c/2 Summation of Energy, TS =58 fps, Plane = c/2 



0 100 200 300 400 0 100 200 300 400 


POD Mode POD Mode 


Summation of Energy, TS =44 fps. Plane = 2c/3 Summation of Energy, TS =58 fps, Plane = 2c/3 



0 100 200 300 400 0 1 00 200 300 400 


POD Mode POD Mode 

Fig. 6 Sumation of Energy per POD mode without 
mean included. 

come as no surprise that it resembles the mean velocity 
profile. 


Eigenvectors in X direction Eigenvectors in Y direction 



y mm y mm 



Fig. 7 Eigenvector Mode 1 at 44/p.s, 15° AOA, c/2 

Reconstruction of a Selected Sample 

With the POD expansion coefficients solved for in 
Equation 4, the flow field can be reconstructed by uti- 
lizing Equation 3. By rebuilding the flow field in this 
manner a minimum number of POD modes can be de- 
termined for adequately representing the velocity field 
of that snapshot. 

To properly determine the minimum number of 
POD modes to be used to rebuild the instantaneous 
snapshot of the flow field, original snapshots of the 
wake flow field need to be used for direct comparison. 
Figure 4 shows instantaneous snapshots of the wake 
flow field at 15° angle of attack, with an airspeed of 
44 fps at a downstream location of c/2. 

By using only 20 out of 2400 POD modes, shown in 
Figure 11, an excellent reconstruction of the instanta- 
neous snapshot of the wake flow field results. This can 
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ficient for all other cases examined. 


X direction 


Reconstruct od U Velocity Flow Field 


Reconstructed V Velocity Flow Field 


Reconstructed W Velocity Flow Field 


Reconstructed V and W Velocity Flow Field 


Fig. 8 Eigenvector Mode 2 at 44 fps, 15° AOA, c/2 


Fig. 11 Rebuild of Instantaneous Snapshot with 
20 POD Mode at a = 15° AOA, airspeed = 44fps, 
plane = c/2 


X direction 


Reconstructed U Velocity Flow Field 


Reconstructed V Velocity Flow Field 


Eigenvectors in Z direction 


Plot of Eigenvectors in Y and Z direction 


Reconstructed W Velocity Flow Field 


Reconstructed V and W Velocity Flow Field 


Fig. 9 Eigenvector Mode 3 at 44 fps, 15° AOA, c/2 


Eigenvectors in X direction 


Fig. 12 Reconstruction of the Instantaneous Snap- 
shot with 2400 POD Modes at a = 15° AOA, air- 
speed = 44 fps, plane = c/2 


Results of Estimation Procedure 

If one considers the distance between the trailing 
edge of the MAV wing to the PIV sampling planes, 
the simultaneous sampling of the strain gauges and 
PIV system need to be phased aligned. Because the 
PIV sampling planes are c/2 and 2c/3 positions down- 
stream from the trailing edge, the flow that is being 
sampled by the PIV in the wake has already passed 
over the wing’s surface. Or in other words, the ed- 
dies that are shed from the MAV wing’s surface have 
to travel a distance downstream before they encounter 
the PIV sampling plane. To phase align these mea- 
surement we define a convective eddy time based on 
Equation 9 


Z direction 


Fig. 10 Eigenvector Mode 4 at 44 fps, 15° AOA 


be seen by direct comparison to the results presented 
in Figure 12 where all 2400 POD modes are included 
in the reconstruction. Twenty modes were found suf- 
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Eigenvectors in Z direction 









Tce = x/U c (9) 

where x is the distance from the middle of the MAY 
wing to the plane and U c is the convective velocity 
that the eddies move at. 20 The convective eddy time 
is estimated to be O.OOSsec for the first plane c/2 where 
U c = 0.7t/oo and Uoo = 44/ps. 

Figure 13 presents a set of simultaneously sampled 
strain gauge measurements along with the PIV TTL 
pulse signal which allows for phase alignment of the 
strain gauge and PIV data. 




Fig. 13 Simultaneously Sampled Strain Gauge and 
TTL Pulse Signals. 


Figure 14 presents an estimate of the wake flow field 
using the estimation procedure described by Equations 
5, 6 and 7 with 20 POD modes. From the results 
presented in Figure 14, the estimate of the wake flow 
field shows that the technique captures much of the low 
dimensional features of the flow. This can be clearly 
seen by comparing these results to those presented in 
Figure 11. 

The focus of this work is on the sensing aspect of the 
control problem. Our proposal for using surface strain 
measurements to estimate the POD expansion coeffi- 
cients in the wake, which are then used to estimate 
the velocity wake flow field when combined with the 
POD eigenfunctions is working to capture the large 
scale features of the flow. 


7 OF 


The POD itself can be viewed as a filter, the esti- 
mation procedure here acts like a filter as well. This 
can be seen by examining the results shown in Fig- 
ure 15 where we include in our estimation all 2400 
POD modes. As can be seen from the figure the es- 
timation of the velocity does not improve with the 
additional POD Modes. In fact comparing the results 
from Figure 15 with those of Figure 11 where only 20 
POD modes were used you can directly see the filter- 
ing effect, of the estimation procedure. The inability of 
the estimation procedure to capture the finer scales is 
probably do to the fact the strain gauges do not pick 
up the signature of the small scales. For control of the 
large scales in the flow however, this loss of informa- 
tion is not critical. 



Fig. 14 Estimation of the Wake Flow Field with 
20 Modes at a = 15°, airspeed = 44 fps, plane = c/2 




40 20 0 


y mm 


LSE Reconstruced V and W Velocity Flow Field 



50 40 30 20 10 0 


y mm 


Fig. 15 Estimation of the Wake Flow Field with 
2400 Modes at a = 15°, airspeed = 44 fps, plane = c/2 
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Conclusions 

A low dimensional tool for flow-structure interaction 
problems based on Proper Orthogonal Decomposition 
(POD) and modified Linear Stochastic Estimation 
(mLSE) was proposed and applied to a Micro Air Ve- 
hicle ( MA V ) wing. The method utilized the dynamic 
strain measurements from the wing to estimate the 
POD expansion coefficients from which an estimation 
of the velocity in the wake can be obtained. 

With the statistical data base that was required for 
the low dimensional tool, the statistical quantities for 
the mean flow field show that the wing tip vortex dom- 
inates the flow. From the POD eigenfunctions, the I ,f 
mode contains 95% or more of the energy in the wake 
flow field. This captures nicely the wing tip vortex. 
With 20 out of 2400 POD modes, a reasonable recon- 
struction of the instantaneous snapshot of the wake 
flow field can be obtained. 

Utilizing the simultaneously sampled strain gauges 
and the wake flow field measurements in conjunction 
with mLSE, an estimation of the wake flow field with 
the lower POD modes, provides a reasonable estima- 
tion of the wake flow field. 

From the results shown, using a low dimensional tool 
for flow-structure interaction provides a valid method- 
ology for estimating the wake flow field from the 
dynamic strain gauges alone. This provides a start- 
ing point for developing control strategies for MAV’s 
where the vehicle itself can morph and essentially be- 
come a part of the control system. 
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